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Abstract — This work presents set-membership adaptive 
algorithms based on time-varying error bounds for CDMA 
interference suppression. We introduce a modified family of 
set-membership adaptive algorithms for parameter estimation 
with time-varying error bounds. The algorithms considered 
include modified versions of the set-membership normalized 
least mean squares (SM-NLMS), the affine projection (SM-AP) 
and the bounding ellipsoidal adaptive constrained (BEACON) 
recursive least-squares technique. The important issue of error 
bound specification is addressed in a new framework that takes 
into account parameter estimation dependency, multi-access and 
inter-symbol interference for DS-CDMA communications. An 
algorithm for tracking and estimating the interference power 
is proposed and analyzed. This algorithm is then incorporated 
into the proposed time-varying error bound mechanisms. 
Computer simulations show that the proposed algorithms are 
capable of outperforming previously reported techniques with a 
significantly lower number of parameter updates and a reduced 
risk of overbounding or underbounding. 

Index Terms — Set-membership filtering, adaptive filters, DS- 
CDMA, interference suppression. 



I. Introduction 

Set-membership filtering (SMF) fTJ, ||2], ||3], ||4l represents 
a class of recursive estimation algorithms that, on the basis of 
a pre-determined error bound, seeks a set of parameters that 
yield bounded filter output errors. These algorithms have been 
used in a variety of applications such as adaptive equalization 
111 and multi-access interference suppression B, id- The 
SMF algorithms are able to combat conflicting requirements 
such as fast convergence and low misadjustment by introduc- 
ing a modification on the objective function. In addition, these 
algorithms exhibit reduced complexity due to data-selective 
updates, which involve two steps: a) information evaluation 
and b) update of parameter estimates. If the filter update 
does not occur frequently and the information evaluation 
does not involve much computational complexity, the overall 
complexity can be significantly reduced. 

The adaptive SMF algorithms usually achieve good con- 
vergence and tracking performance due to an adaptive step 
size for each update, and reduced complexity resulting from 
data selective updating. However, the performance of SMF 
techniques depends on the error-bound specification, which 
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is very difficult to obtain in practice due to the lack of 
knowledge of the environment and its dynamics. In wireless 
networks characterized by non-stationary environments, where 
users often enter and exit the system, it is very difficult to 
choose an error bound and the risk of overbounding (when the 
error bound is larger than the actual one) and underbounding 
(when the error bound is smaller than the actual one) is 
significantly increased, leading to performance degradation. 
In addition, when the measured noise in the system is time- 
varying and the multiple access interference (MAI) and the 
intersymbol interference (ISI) encountered by a receiver in 
a communication system are highly dynamic, the selection 
of an error-bound is further complicated. This is especially 
relevant for low-complexity estimation problems encountered 
in applications such as mobile terminals and wireless sensor 
networks [8], where the sensors have limited signal processing 
capabilities and power consumption is of central importance. 
These problems suggest the deployment of mechanisms to 
automatically adjust the error bound in order to guarantee 
good performance and a small update rate (UR). It should also 
be remarked that most of prior work on adaptive algorithms 
for interference suppression llTSI is restricted to systems with 
short codes. However, the proposed adaptive techniques are 
also applicable to systems with long codes provided some 
modifications are carried out. For downUnk scenarios, the 
designer can resort to chip equalization ^241 followed by a 
despreader For an uplink solution, channel estimation algo- 
rithms for aperiodic sequences ll25l . ||26] are required and the 
sample average approach for estimating the covariance matrix 
R = E[r{i)r^ (i)] of the observed data r(i) has to be replaced 
by R = PP^ + (T^I, which is constructed with a matrix P 
containing the effective signature sequence of users and the 
variance cr^ of the receiver's noise ||27|. 

Previous works on time-varying error bounds include the 
tuning of noise bounds in |9|, [lOJ, the approach in f\3\ 
which assumes that the "true" error bound is constant, and 
the parameter-dependent error bound recently proposed in 
[|1 U . [12J with frequency-domain estimation algorithms. The 
techniques so far reported do not introduce any mechanism for 
tracking the MAI and the ISI and incorporating their power es- 
timates in the error bound. In addition, the existing approaches 
with time-varying bounds have not been considered for more 
sophisticated adaptive filtering algorithms such as the affine 
projection (AP) and the least-squares (LS) based techniques. 

In this work, we propose and analyze a low-complexity 
framework for tracking parameter evolution and MAI and ISI 
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power levels, that relies on simple channel and interference 
estimation techniques, and encompasses a family of set- 
membership algorithms [2], [3), [6], [141 with time- varying 
error bounds. Specifically, we present modified versions of the 
set-membership normalized least mean squares (SM-NLMS) 
m, the affine projection [31 (SM-AP) and the bounding 
ellipsoidal adaptive constrained (BEACON) JS], lfT4l recur- 
sive least-squares (RLS) algorithm for parameter estimation. 
Then, we incorporate the proposed mechanisms of interference 
estimation and tracking into the time-varying error bounds. 
In order to evaluate the proposed algorithms, we consider a 
DS-CDMA interference suppression application and adaptive 
linear multiuser receivers in situations of practical interest. 

This work is organized as follows. Section II briefly de- 
scribes the DS-CDMA system model and linear receivers. 
Section III reviews the SMF concept with time-varying error 
bounds and is devoted to the derivation of adaptive algorithms. 
Section IV presents the framework for time-varying error 
bounds and the proposed algorithms for channel, interference 
estimation and tracking. Section V is dedicated to the analysis 
of the algorithms for channel, amplitude, interference estima- 
tion and their tracking. Section VI shows and discusses the 
simulations results, while Section VII gives the conclusions. 

II. DS-CDMA System Model and Linear Receivers 

Let us consider the downlink of a symbol synchronous DS- 
CDMA system with K users, N chips per symbol and Lp 
propagation paths 1 18 1. We assume that the delay is a multiple 
of the chip rate, the channel is constant during each symbol 
interval and the spreading codes are repeated from symbol 
to symbol. The received signal r{t) after filtering by a chip- 
pulse matched filter and sampled at chip rate yields the M- 
dimensional received vector 
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where M = N + Lp — 1, n[i] = [ni[i] . . . nM[i\]'^ is the 
complex Gaussian noise vector with zero mean and covariance 
matrix i?[n[i]n^[i]] = cr^I, where (•)^ and (•)^ denote 
transpose and Hermitian transpose, respectively. The quantity 
E[-] stands for expected value, the user k symbol is bk[i] 
and is assumed to be drawn from a general constellation. 
The amplitude of user k is Ak[i] and ri[i] is the intersymbol 
interference (ISI) for user k. The M x Lp convolution matrix 
C/c that contains one-chip shifted versions of the signature 
sequence for user k expressed by = [ak{l) ■ ■ ■ ak{N)]'^ 
and the Lp x 1 vector h[i] with the multipath components are 
described by: 
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In this model, the ISI span and contribution rii^[i] are functions 
of the processing gain N and Lp. If 1 < Lp < N then 3 



symbols would interfere in total, the current one, the previous 
and the successive symbols. In the case of TV < Lp < 2N then 
5 symbols would interfere, the current one, the 2 previous and 
the 2 successive ones. In most practical CDMA systems, we 
have that 1 < Lp < iV and then only 3 symbols are usually 
affected. The reader is referred to UMTS channel models [22], 
which reveal that the channel usually affects at most 3 symbols 
(it typically spans a few chips). 

The multiuser linear receiver design corresponds to deter- 
mining an FIR filter Wk[i] — [wa[i] wi[i] . . . wm-i] with 
M coefficients that provides an estimate of the desired symbol 
as given by 



bkii] = sg 
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where the quantity 5R(-) selects the real part and sgn(-) is 
the signum function. The quantity Zk[i] = wj^[i]r[i] is the 
output of the receiver parameter vector wj. for user k, which 
is optimized according to a chosen criterion. 

III. Set-Membership Adaptive Algorithms with 
Time- Varying Error Bounds and Problem 
Statement 

In this section, we describe a framework that encompasses 
modified set-membership (SM) adaptive algorithms with time- 
varying error bounds for communications applications. In an 
SM filtering |2| framework, the parameter vector w/j[«] for 
user fc in a multi-access system is designed to achieve a 
specified bound on the magnitude of the estimation error 
efc[i] = bk[i] — w^[i]r[i]. As a result of this constraint, the 
SM adaptive algorithm will only perform filter updates for 
certain data. Let represent the set containing all Wk[i] 

that yields an estimation error upper bounded in magnitude by 
a time-varying error bound 7^ [i] . Thus, we can write 

©fcW = n {^fe ^ -I 1^ '^^[*]}' (4) 

where r[i] is the observation vector, Sk is the set of all possible 
data pairs {bk[i], r[i]) and the set 9fc[«] is referred to as 
the feasibility set for user k, and any point in it is a valid 
estimate Zk[i] = wj^[i]r[i]. Since it is not practical to predict 
all data pairs, adaptive methods work with the membership 
sets ^pk^i — nL=i ^fc,ni provided by the observations, where 
'Hk.m = {wfe e C*^ : \bk\m\ - Zfc[m]| < 7fc[m]} is the 
constraint set. It can be seen that the feasibility set 
is a subset of the exact membership set at any given time 
instant. The feasibility set is also the limiting set of 

the exact membership set, i.e., the two sets will be equal 
if the training signal traverses all signal pairs belonging to 
Sk- The idea of the SM algorithms is to adaptively find 
an estimate that belongs to the feasibility set 0fc[i]. One 
alternative is to apply one of the many OBE algorithms such as 
the bounding ellipsoidal adaptive constrained (BEACON) 1 14|, 
|l6l recursive least-squares (RLS) algorithm, which tries to 
approximate the exact membership set with ellipsoids. Another 
way is to compute a point estimate through projections using, 
for example, the information provided by the constraint set 
T-ik,!, as done by the set-membership NLMS (SM-NLMS) 
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and the affine projection |]3] (SM-AP) algorithms. In order 
to devise an effective SM algorithm, the error bound •jki'i] 
must be appropriately chosen. Due to time-varying nature of 
many practical environments, this error bound should also 
be adaptive and adjustable to certain characteristics of the 
environment for the SM estimation technique. The natural 
question that arises is: how to design an efficient and effective 
mechanism to adjust 7fe[«]? In what follows, we will present a 
modified family of SM adaptive algorithms that rely on general 
time-varying error bounds. Specifically, we will consider the 
SM-NLMS SM-AP fj] and BEACON (14), f6l algorithms 
and we will modify them such that they will operate with 
general time-varying error bounds. 

A. SM-NLMS Algorithm with Time-Varying Bounds 

In order to derive an SM-NLMS adaptive algorithm with 
time-varying bounds using point estimates, we consider the 
following optimization problem 

minimize ||wfe[i + 1] — Wfe[z]|p 

subject to {bk[i\ - wf [i + l]r[i]) = gk\i] 

In order to solve the above constrained optimization problem, 
we resort to the method of Lagrange multipliers ||T], Il23l , 
which yields the unconstrained cost function 



The SM-AP adaptive algorithm with time-varying bounds can 
be derived from the optimization problem 

|2 

(12) 



minimize ||wi.[i + 1] — Wfe[i] 
subject to {hk[i\ - Y"\i]wk[i + 1]) = gk[i\- 



In order to solve the above problem, we employ the method 
of Lagrange multipliers and consider the unconstrained cost 
function 



= ||wfc[z + 1] - yvk[i] 
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(bfeH-Y«Hwfe[i + l]-g,H)^A 
(13) 



C = ||wfc[i + 1] - yvk[i] 
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where * denotes complex conjugate, A is a Lagrange multiplier 
and gk[i] is the time-varying set-membership constraint for 
user k. Taking the gradient terms of (6) with respect to Wk[i + 
1] and A*, and setting them to zero, leads us to a system of 
equations. Solving these equations yields: 



where A is the vector with Lagrange multipliers and gfe[i] 
is a constraint vector By calculating the gradient terms of 
(13) with respect to Wk[i + 1] and A, setting them to zero 
and solving the resulting equations we arrive at the following 
algorithm: 

tfeW - (Y^HYW + <5I)-i(e,H - gW), (14) 

Wfe[i + 1] =Wfe[z]+YHtfc[z], (15) 

where (5 is a small constant inserted in addition to the term 
Y-^[i]Y[i] for improving robustness. If we select ek[i] — 
gk\i] = (efc[?] - ■^sgn{ek[i\))u = (1 - 7fe[«]/|efe[i]|)efc[i]u, 
where the a posteriori errors [i — j] are kept constant for 
1, . . . , P - 1 and u = [1 . . . 0]^, we obtain the 
6wing recursion for the update of tk[i]'- 

tk[i] = (Y^WYH + SI)-\1 - 7/|efe[z]|)efeWu. (16) 



Gk [l_ 



= bk H] - 



w.. i r i 



Wfc[z + 1] - WfeW + (r^[i]rH)-i(efcH - 5feW)*rH, 



(7) 

(8) 



Substituting (16) into (15) and using the bound constraint, we 
obtain the following SM-AP algorithm: 

i]Y[2]{Y"[2]Y[i]+SI)-^ek[i]u, (17) 

(18) 



(1 -7fcW/|eA;K 




if \ek\i]\ > 7fcM 
otherwise. 



where ek[i] is the error for user k. By choosing gk[i] such that 
efc[z] lies on the closest boundary of 0fc[z] and considering a 
time-varying error bound 7fc[i], i.e., gk[i] = 'yk[i\sgn{ek[i]) 
Q, we obtain the following data dependent update rule and 
step size 



^k[i - 
1 



— J r 



mi] 



1] = Wk[i] + Hw[i]el[i]r[i 
(l-7fcW/|e:.M 



(9) 



otherwise. 



(10) 



B. SM-AP Algorithm with Time-Varying Bounds 

In order to describe a modified SM-AP algorithm with 
time-varying bounds, let us first define the observation matrix 
Y[i] — [r[i] ... r[i — P + 1]], the desired output vector 
bfe[i] = [bk[i] ■ ■ ■ bk[i — P + 1]]"^ that comprises P outputs 
and the error vector 



ekM 



-r»[z-P+l]wfc[z] 



-Y^ 



l\Wk[l\ 



(11) 



The SM-AP algorithm described here has computational com- 
plexity of UR X 0{PM + 2KinvP'^), where Ki^v is a factor 
required to invert a PxP matrix (T\ and UR is the update rate. 
Note that the SM-AP is a generalized case of the SM-NLMS 
where P data vectors are used to increase the convergence 
speed. 

C. BEACON Adaptive Algorithm with Time-Varying Bounds 

Here, we propose a modification for a computationally 
efficient version of an optimal bounding ellipsoidal (OBE) al- 
gorithm called the Bounding Ellipsoidal Adaptive Constrained 
Least-Squares (BEACON) algorithm [6|, which is closely 
related to a constrained least-squares optimization problem. 
The proposed technique amounts to deriving the BEACON al- 
gorithm equipped with time-varying bounds. Unlike the other 
previously reported low-complexity algorithms ||2], ||3] and the 
modified SM-NLMS and SM-AP techniques described in the 
previous subsections, the modified BEACON recursion has 
the potential to achieve a very fast convergence performance, 
which is relatively independent from the eigenvalue spread 
of the data covariance matrix as compared to the stochastic 
gradient algorithms [jT]. 



4 



The proposed BEACON algorithm with a time-varying 
bound can be derived from the following optimization problem 

m 



minimize ^ A' - w^[i]r[i 



1=1 



(19) 



subject to \hk[i] — w k[i]^ r[i]\'^ = 



Wfe i 



Wfe i 



1 + Xk\i]tk\i] 

where the prediction error is = bk[i] — w^[i — l]r[«], 

tf,[i] = r^mPkii - and Xk[i] > 0. In order to 

compute the optimal value for Xk[i], the algorithm considers 
the following cost function ||6l: 



C 



Ml 



1 + Xk\i]tk[i] 
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The minimization of the cost function in (24) leads t 
innovation check of the proposed BEACON algorithm: 



XM = 




1 



if iaMI>7feM 
otherwise. 



A. Parameter Dependent Bound 

Here, we describe a parameter dependent bound (PDB), 
that is similar to the one proposed in [TTl and considers the 
evolution of the parameter vector Wfe[z] for the desired user 
(user k). The proposed PDB recursion computes a bound for 
SM adaptive algorithms and is described by: 



The constrained problem above can be recast as an uncon- 
strained one and be solved via an unconstrained least-squares 
cost function using the method of Lagrange multipliers given 
by 



7fc[i + !] = (!- /3)7fe[z] + P^Ja\\vfk[^\\\ 



(26) 



(20) 

where Afe [/] plays the role of Lagrange multiplier and forget- 
ting factor at the same time for user k. The solution to the 
above optimization problem is obtained by taking the gradient 
terms with respect to w/j[«] and making them equal to zero. 
After some mathematical manipulations we have 

WfcH =Wfe[i~l] + Afe[z]PfcWrH(6feH-wf [z]r[z])*. (21) 

By using the constraint hk[i] — 'w^[i]v[i] — 

(CfcM/l4WI)7feW = aw - (A.MGfc[z](4W/ie;cWI)7fcW), 

and the matrix inversion lemma [JJ, we can arrive at the 
BEACON algorithm with time-varying bounds described by 

A,WP,[z-l]rWr^HPfe[i-l]^ ^22) 

(23) 



(24) 



IV. Algorithms for Time- Varying Error Bour 
Interference Estimation and Tracking 

This section is devoted to time-varying error bounds th 
corporate parameter and interference dependency. We pn 
a low-complexity framework for time-varying error bo 
interference estimation and tracking. A simple and eff 
algorithm for taking into consideration parameter depen( 
is introduced and incorporated into the error bound. A proce- 
dure for estimating MAI and ISI power levels is also presented 
and employed in the adaptive error bound for SM algorithms. 
The proposed algorithms are based on the use of simple rules 
and parameters that behave as forgetting factors, regulate the 
pace of time averages and selectively weigh some quantities. 



where /3 is a forgetting factor that should be adjusted to ensure 
an appropriate time-averaged estimate of the evolutions of the 
parameter vector Wfc[z], Q;||wfc[i]|p(7y[z] is the variance of the 
inner product of Wfc[i] with n[z] that provides information on 
the7ifi|))ution of Wfe[i], a is a tuning parameter and (7^[i] is 
an estimate of the noise power. This kind of recursion helps 
avoiding too high or low values of the squared norm of w^: [i] 
and provides a smoother evolution of its trajectory for use in 
the time-varying bound. The noise power at the receiver should 
be estimated via a time average recursion. In this work, we 
will assume that it is known at the receiver 



B. Parameter and Interference Dependent Bound 

In this part, we develop an interference estimation and 
tracking procedure to be combined with a parameter dependent 
bound and incorporated into a time-varying error bound for 
SM recursions. The MAI and ISI power estimation scheme, 
outlined in Fig. 2, employs both the RAKE receiver and the 
hnear receiver described in (3) for subtracting the desired user 
signal from v[i] and estimating MAI and ISI power levels. 
With the aid of adaptive algorithms, we design the linear 
receiver, estimate the channel modeled as an FIR filter for the 
RAKE receiver and obtain the detected symbol bk[i], which is 
combined with an amplitude estimate Ak [i] for subtracting the 
desired signal from the output Xk[i\ of the RAKE. Then, the 
difference dk[{\ between the desired signal and Xk[i\ is used 
to estimate MAI and ISI power. 
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Fig. 1. Block diagram of the proposed interference estimation and tracking 
algoritiim. 
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Channel Estimation: 

Let us first present a simple channel estimation algorithm for 
designing the RAKE receiver. Consider the constraint matrix 
Cfc that contains one-chip shifted versions of the signature 
sequence for user k defined in (2) and the assumption that the 
symbols bk \i] are independent and identically distributed (i.i.d) 
and statistically independent from the symbols of the other 
users. The proposed channel estimation algorithm is based on 
the following cost function 



C{h[{\,Ak[{\) = E[\\AkW>k[i]Ckh[i 

= E[\\AkW>kWk[^-v[i 



m\ 

21 



(27) 



where h[i] is an estimate of the channel h[i\ and ffc[i] = 
Cfch[i] is the RAKE receiver for user k with the estimated 
channel. By taking the gradient terms of (27), making them 
equal to zero, we can devise a stochastic gradient (SG) channel 
estimation algorithm as follows: 

h[i + 1] = h[i] - ^XhAk[l\C^bl\i]{bk[{\Ck^i\ - v\i]). (28) 

AmpUtude Estimation: 

The amplitude has also to be estimated at the receiver in 
order to provide this information for different tasks such as 
interference cancellation and power control. The proposed in- 
terference estimation and tracking algorithm needs some form 
of amplitude estimation in order to proceed with the estimation 
of the interference power To estimate the amplitudes of the 
associated user signals, we describe the following optimization 
problem with the cost function in (27) 



Ak[i] = argminC(h[i], 



arff min E\ 



\Ak[i]bk[i\Ckh.[i] - v\i\ 



(29) 



where f^[z]sfe[i] = pk[i] and [i]sj[i] = pi,j[i] for j ^ 
1. The symbol pk represents the cross-correlation (or iimer 
product) between the effective signature and the RAKE with 
perfect channel estimates. The symbol represents the 

cross-correlation between the RAKE receiver and the effective 
signature of user j. The second-order statistics of the output 
of the RAKE in (31) are described by: 



E[\xk[ 



Al\i]pl^]E[\bM' 
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H~ -Hi 



1=1 



+ f^E[r,\^]r,"mk + ^^E[n\^]n"[ 



(32) 

From the analysis above, we conclude that through the 
second-order statistics one can identify the sum of the power 
levels of MAI, ISI and the noise terms. Therefore, our strategy 
is to obtain instantaneous estimates of the MAI, the ISl and 
the noise from the output of a RAKE receiver, subtract the 
detected symbol in (3) from this output (using the more reh- 
ab le multiuser receiver (wfe[i])) and to track the interference 
(MAI + ISI + noise) power as shown in Fig. 2. Let us define 
the difference between the output of the RAKE receiver and 
the detected symbol for user 1: 

K 

dk[i] =Xk[i] -Ak[i]bk[i] ^ Ak[i]bk[i\ik [i]sk 

k=2 




MAI 



i\n\i\ 



In order to solve the problem above efficiently we describe a 
simple SG algorithm to estimate the amplitude of user k, as 
given by 

+ =ifc[i] -^^6^[i]h^[i]Cf (ife[^]6fe[^]Cfch[^] -r[i]). 



(33) 



By taking expectations on and taking into account the 

assumption that MAI, ISI and noise are uncorrelated we have: 



(30) E[\dk[i]\' 



Interference Estimation and Tracking: 



K 

k=2 



k [i% 



(34) 



+ i^[i\E[ri\i\ri"mm+<^'i{lWk[ 



Let us consider the RAKE receiver with perfect channel 
knowledge, whose parameter vector fk\i] = Ckh[i] for user k 
(desired one) estimates the effective signature sequence at the 
receiver, i.e. Cfeh[i] = Sk[i]. The output of the RAKE receiver 
is given by: 



where the above equation represents the interference power. 
Based on time averages of the instantaneous values of the 
interference power, we introduce the following algorithm to 
estimate and track ii^[|dfe[i]p]: 



Xk\l 



= Ak[i\bk[i]fk[i\sk\ 



i + 



desired signal 



K 

i=2 



v[i + i] = ii-p)v\i] + p\dk[i\\'' 



(35) 



i]bj[i\tf 



i\s. 



MAI 



ISI 



(31) 



[*Jwhere /3 is a forgetting factor. To incorporate parameter 
dependency and interference power for computing a more 

—effective bound, we propose the parameter and interference 
dependent bound (PIDB): 

7fe[i + 1] = (1 - Phk\i] +/3(Vr v^\i] + ^/a\\wkK[i]), 

(36) 
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where v[i] is the estimated interference power in the multiuser 
system and r is a weighting parameter that must be set. 
Similarly to (26), the equations in (35) and (36) are time- 
averaged recursions that are aimed at tracking the quantities 

|c?fe[i]p and {\/Tv'^[i\ + \J a||wfc||2(f„^[i]), respectively. The 
equations in (35) and (36) also avoid undesirable too high 
or low instantaneous values which may lead to inappropriate 
time-varying bound 7^ [i] . 

V. Analysis of The Algorithms 

In this section we analyze the channel estimation and 
interference estimation algorithms described in the previous 
section. We focus on the convergence properties of the al- 
gorithms in terms of the step size parameters /i/j, and /3 
used for the channel, the amplitude and the interference power 
estimators, respectively. 

A. Channel Estimator 

In order to analyze the SG channel estimator given in (28), 



algorithm, the magnitude of this geometric ratio must be less 
than 1 for all fc (|1 — ^/iCr^^f7^C|^Cfc| < 1). This means that 
-. The reviewer is referred to [23] 



let us first define an error vector e/, [i] = h[i] — hopt- 
subtracting hopt from the SG recursion in (28), we get 

e,J^ + 1] = e„W - M^ifc^Cf &t4*](A-WfefeWCfehH 



By 



-r[z]) 



opt. 



-r[z]) 

= [1- ^ihAk[l]\hkWC^Ck\eh[i] 

- fihAk \i]bl [i\ Cf (ife [i\bk \i]Ckhopt 
= [l-^„ifc[z]|6fe[z]pCfCfe]e;,H 



(37) 



By considering that the error vector eh[i], ek.opt[i], the sig- 
nal components X^feLi-^^-^W from the data vector r[i] = 
J2^=i^k[i] + n[i] given by (1) and n[i] are statistically 
independent and computing the covariance matrix of the error 
vector, i.e. K[z] = ii^[e/i[i]e^[i]], we obtain 

K[i + 1] = [l-fi^aiy,C^Ck]K[i][l-f,Haiy,C^Ck] 

(38) 

where = E[\bk\i]\% ^1 - E[\Ak\i]\^] and MSE„,„ = 
-E'[||efc,opt[i]|P] stands for the minimum MSE achieved by the 
estimator The recursive rule/algorithm in (28) is asymptot- 
ically unbiased and will converge to the optimum channel 
estimator hopt if the step size fih satisfies the following 
condition ^ 

0<^ih< ^ , (39) 



where Xmax is the largest eigenvalue of the matrix C^Cfc. 
The condition above with concern to the step size /i/j arises 
from difference equations. The quantities generated in (38) 
represent a geometric series with a geometric ratio equal to 
(1 — iiha\ a^C^ Ck)- For stability or convergence of this 



< A^'' < 

k 

for further details. 



B. Amplitude Estimator 

In order to analyze the SG amplitude estimator described in 
(30), let us first define an error signal eA,k[i] = ^fcH —Ak^opt- 



By subtracting Af 



<^A, [j + 1] = <^A, [i 

^ e^, [i 
+ Ak. 
= <^A, [i 



opt from the equation in (30) we obtain 



- MA^^Wh^'WCf (ifc[z]6fc[i]Cfeh[z] - r[i]) 

~f,Abl[l]h"[l]Ck{i^A,kll] 

opt)bk[i]Ckh[{\ - r[i]) 
-MA|&fcWI'hHCfCfch[z]eA,W 
- fiAbmh"[i]C^ {Ak,optbk[^Ckh[i] - r[z]) 
= (l-MA|6feMphMCfCfeh[z])e^,H 

(40) 

By considering that the error e^W, ©yi.optW = 
Ak,optbk[i\Ckh[{\ - r[i\, the signal components Y.k=i^k[i\ 
from the data vector r[i] = X^fLi + '^M given by (1) and 
n[i] are statistically independent, and by computing the mean- 
squai-ed error, i.e., MSE(Afc[i])[i] = i^Aji] = i?[|e[i]p], we 
get 

, 2 , 



MSE(AfeH)[i + l] 



(l-MACT,'i?[||hf[z]Cf|p])^X^Jz] 

S[||hfHCf ||2]MSE^,™„, 

(41) 



2 2 



where MSE^.m^n = eA,opt[*]|| ]■ The cross multiplication 
between the terms will vanish as a result of the statistical 
independence between them. The general algorithm in (30) 
will converge asymptotically and in an unbiased way to the 
Ak,opt provided the step size fiA is chosen such that 

'^'^^mmmfY '''' 

The above range of values has to be tuned in order to ensure 
good convergence and tracking of the amplitude. 

C. Interference Estimation and Tracking Algorithm 

Let us describe in a general form the mechanisms for 
interference estimation and tracking given in (26) and (36). 
We can write without loss of generality 



7fe[i + l] = (l-/3)7fcW + /3PoM 



(43) 



where Po[i] can account either for a parameter dependent 
bound (PDB), as described in Section IV. A, or for a parameter 
and interference dependent bound (PIDB), as the one detailed 
in Section IV.B. 

Our goal is to establish the convergence of a general 
stochastic recursion, as the one given in (43), in terms of the 
mean-squared error (MSE) at iteration [i], as described by 

|21 



MSE(7[z])W = i?[|7,[ 



7fc,opt| 



(44) 
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where 7fe,opt is the optimal parameter estimate for 7fc[i]. 

We will show that under certain conditions on f3, the 
sequence converges to the optimal 7fe^opt in the mean-square 
sense. Let us define the error 

e'f[i\ = 7fcW - 7fe,opt, 
and substitute the above equation into (44), which yields 

[z + 1] = (1 - /3)e^ W + /3(PoM - 7/c.op,) 
-(l-/3)e^H+/3efe,optW, 
The MSE at time instant [i + 1] is given by 

MSE(7[i + + l]=e;[i + l]e^[i + 1] = \e^[i + l]p 

+ /3^ efc,opt[i]efc^opt[i] 
+ {1-/3)13 e^WeLptW 
+ {1-/3)13 e;[z]efc,op,[z]. 

(46) 

By taking the expected value on both sides and assuming that 
e-y[i] and Cfc^optW are statistically independent we have 

E[MSE[i + l]] = E[\e-y[i + 1]\^] 
= Ky[i + l] 

= (l-/3)if^W(l-/3) + /32lefc,„pt[z]|2, 

(47) 

where |efc optW|^ is the minimum MSE achieved by the esti- 
mator 

We can notice that the cross-multiplication terms will vanish 
as a result of the statistical independence between the terms. 
The general recursive rule in (43) will converge asymptotically 
provided the step size /3 is chosen such that 

< /3 < 2. (48) 

The above range of values has to be adjusted in order to ensure 
good convergence and tracking of the parameter dependency 
and/or the interference modeled here as the quantity Po [i] . 

VI. Simulations 

In this section we assess the performance of the proposed 
and existing adaptive algorithms in several scenarios of prac- 
tical interest: 

. The NLMS, the SM-NLMS 0, the SM-NLMS of Guo 
and Huang iTTl and the proposed SM-NLMS with the 
parameter-dependent (PDB) and parameter and interfer- 
ence dependent (PIDB) time- varying bounds. 
. The AP |1|, the SM-AP |3| and the proposed SM-AP 

with the PDB and PIDB time-varying bounds. 
. The BEACON M and the proposed BEACON algorithms 

with the PDB and PIDB time-varying bounds. 
We consider for the sake of simplicity binary phase-shift- 
keying (BPSK) modulation, a DS-CDMA system with Gold 
sequences of length = 31 and typical fading channels 
found in mobile communications systems that can be modeled 
according to Clarke's model 121]. The channels experienced 
by different users are identical since we focus on a downlink 



scenario and the desired receiver processes the transmissions 
intended to other users (MAI) over the same channel as its own 
signal. The channel coefficients are hi [i] — piai [i], where ai [i] 
(I — 0,1, Lp — 1) are obtained with Clarke's model [21]. 
In particular, we employ standard SG adaptive algorithms for 
channel estimation and RAKE design in order to concentrate 
on the comparison between the analyzed algorithms for re- 
ceiver parameter estimation. We show the results in terms of 
the normalized Doppler frequency fdT (cycles/symbol) and 
use three-path channels with relative powers given by 0, —3 
and —6 dB, where in each run the spacing between paths is ob- 
tained from a discrete uniform random variable between 1 and 
2 chips. The channel coefficients hi[i] (I = 0,1, Lp — 1) 
are constant during each symbol interval and change according 
to Clarke's model over time. Since the channel is modelled as 
an FIR filter, we employ a channel estimation filter with 6 
taps as an upper bound for the experiments. Note that the 
delays of the channel taps are multiples of the chip rate and 
are random. Their range coincides with the maximum length 
of the estimation filter which is 6 taps. 

The parameters of the algorithms are optimized and shown 
for each example, the system has a power distribution for 
the interferers that follows a log-normal distribution with 
associated standard deviation of 3 dB and experiments are 
averaged over 100 independent runs. The receivers are trained 
with 200 symbols and then switch to decision-directed mode 
for each data packet. We address the dynamic channel by 
adjusting the receiver weights with the training sequences 
(with length equal to 200 symbols) and then we exploit the 
decision-directed mode to track the channel variations. If the 
channel does not vary too fast then the adaptive receivers can 
track it with this scheme, as will be shown in what follows. 

A. Interference Estimation and Tracking 

In order to evaluate the effectiveness of the proposed 
interference estimation and tracking algorithms, that are in- 
corporated into the time-varying error bounds for tracking the 
MAI and ISI powers, we carried out an experiment, depicted 
in Fig. 2. In this scenario, the proposed algorithm estimates 
of the MAI and ISI powers are compared to the actual 
interference power levels that are generated by the simulations. 
The results obtained show that the proposed algorithms are 
very effective for estimating and tracking the interference 
power in dynamic environments. Specifically, the estimation 
error does not exceed in average 5% of the estimated power 
level, as depicted in Fig. 2. Indeed, the algorithms are capable 
of accurately estimating the interference power levels and 
tracking them, which can be verified by the proximity between 
the curves obtained with the estimation algorithms and the 
actual values. 

B. SINR Performance 

In this part, the performance of the proposed algorithms is 
assessed in terms of output signal-to-interference-plus-noise 
ratio (SINR), which is defined as 

SINR z (49) 
w" [i\Ri [i\w[i\ 
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where Rs [i] is the autocorrelation matrix of the desired signal 
and Ri[i] is the cross-correlation matrix of the interference 
and noise in the environment. The goal here is to evaluate 
the convergence performance of the proposed algorithms with 
time-varying error bounds for the modified SM-NLMS, SM- 
AP and BEACON techniques. Specifically, we consider exam- 
ples where the adaptive receivers converge to about the same 
level of SINR, which illustrates in a fair way the speed of 
convergence of the proposed algorithms and the existing ones. 
We also measure the update rate (UR) of all the SM-based 
algorithms as an important complexity issue. 

The SINR convergence performance of NLMS, AP and 
BEACON algorithms is illustrated via computer experiments 
in Figs. 3, 4 and 5, respectively. The curves in Fig. 3 
show that the proposed SM-NLMS algorithms with the PIDB 
time-varying error bounds achieve the fastest convergence, 
followed by the proposed SM-NLMS-PDB, the SM-NLMS- 
Huang m, the conventional SM-NLMS |2| and the NLMS 
0] recursions. Even though the proposed SM-NLMS-PIDB 
and SM-NLMS-PDB algorithms enjoy the fastest convergence 
rates, they exhibit remarkably lower UR properties, saving 
significant computational resources and being substantially 
more economical than the conventional SM-NLMS algorithm. 

By observing the results for the AP and the BEACON 
algorithms, shown in Figs. 4 and 5, one can notice that the 
results corroborate those found for the NLMS algorithms. 
It should be noted that despite their higher complexity than 
NLMS algorithms, the AP and BEACON techniques have 
faster convergence, better SINR steady-state performance and 
lower URs. 

C. BER Performance 

In this subsection, we focus on the bit error rate (BER) 
performance of the proposed algorithms. We consider a sim- 
ulation setup where the data packets transmitted have 1500 
symbols and the adaptive receivers and algorithms are trained 
with 200 symbols and then switch to decision-directed mode. 



N=31, K=10 users, E /N =15 dB, f ,T=0.00005 

DO d 







f 






— e — AP(P=3) 




— H — SM-AP(P=3)-UR=24.2% 




— ft — SM-AP(P=3)-PDB-UR=20.5% 




i SM-AP(P=3)-PIDB-UR=15.7% 



500 1000 1500 

Number of received symbols 



Fig. 4. SINR Performance of AP algorithms at Eb/No = 15 dB with 
13 = 0.05, a = 8 and r = 2 . 



in which they continue to adapt and track the channel varia- 
tions. 

Firstly, we consider a study of the BER performance and the 
impact on the UR of the fading rate of the channel (fdT) in 
the experiment shown in Fig. 6. We observe from the curves in 
Fig. 6 (a) that the new BEACON algorithms obtain substantial 
gains in BER performance over the original BEACON [6 1 and 
the RLS algorithm [l] for a wide range of fading rates. In 
addition, as the channel becomes more hostile the performance 
of the analyzed algorithms approaches one another, indicating 
that the adaptive techniques are encountering difficulties in 
dealing with the changing environment and interference. With 
regard to the UR, the curves in Fig. 6 (b) illustrate the 
impact of the fading rate on the UR of the algorithms. 
Indeed, it is again verified that the proposed BEACON-PDB 
and BEACON-PIDB algorithms can obtain significant savings 
in terms of UR, allowing the mobile receiver to share its 
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(a) N=31, K=8 users, f T=0.0001 



(b) N=31, E /N =15 dB, f T=0.0001 




e — RLS 

B — BEACON-UR=9.8% 
* — BEACON-PDB-UR=7.3% 
BEACON-PIDB-UR=5.8% 



500 1000 
Number of received symbols 



1500 




-e — AP (P=3) 
-B — SM-AP (P=3) 

— SM-AP (P=3)-PDB 
SM-AP (P=3)-PIDB 



12 16 20 2 4 6 8 10 12 14 16 18 20 



E,/N„ (dB) 



Number of users (K) 



Fig. 5. SINK Perforaiance of RLS and BEACON algorithms at Et/No = 15 Fig. 7. BER performance versus (a) E^/No dB (b) K for AP algorithms 



dB with B = 0.05, a = 5 and r = 1.5 . 



with P = 3. 



-G — RLS 

-a — BEACON 

BEACON-PDB 

—I — BEACON-PIDB 




10 10" 
f J (cycles/symbol) 



10 10"' 
f T (cycles/symbol) 



Fig. 6. (a) BER performance versus f^T and (b) Update Rate (UR) versus 
fdT for the BEACON algorithms for a forgetting factor A = 0.997 for the 
conventional RLS algorithm. The proposed BEACON-PDB abd BEACON- 
PIDB algorithms use /3 = 0.05, a = 5 and t = 1.5. 



more users as compared to the PDB technique for the same 
BER performance. 

VII. CONCLUSIONS 

We proposed SM adaptive algorithms based on time-varying 
error bounds. Adaptive algorithms for tracking MAI and ISI 
power and taking into account parameter dependency were 
incorporated into the new time-varying error bounds. Simu- 
lations show that the new algorithms outperform previously 
reported techniques and exhibit a reduced number of updates. 
The proposed algorithms can have a significant impact on 
the design of low-complexity receivers for spread spectrum 
systems, as well as for future MIMO systems employing either 
CDMA or OFDM as the multiple access technology. The 
proposed algorithms are especially relevant to future wireless 
cellular, ad hoc and sensor networks, where their potential 
to save computational resources may play a significant role 
given the limited battery resources and processing capabilities 
of mobile units and sensors. 



processing power with other important functions and to save 
battery life. 

The BER performance versus the signal-to-noise ratio 
(Efj/No) and the number of users (K) is illustrated in Figs. 7 
for the AP with P — 3. The results confirm the excellent per- 
formance of the proposed PIDB time-varying error bound for a 
variety of scenarios, algorithms and loads. The PIDB technique 
allows significantly superior performance while reducing the 
UR of the algorithm and saving computations. We can also 
notice that significant performance and capacity gains can be 
obtained by exploiting data reuse. From the curves it can be 
seen that the proposed PIDB mechanism with the SM-AP with 
P = 3 can save up to 3 dB and up to 4 dB as compared to 
the PDB and to the SM-AP with fixed bounds, respectively, 
for the same BER performance. In terms of system capacity, 
we verify that the PIDB approach can accommodate up to 4 
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